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Abstract. In the framework of the Chandrasekhar-mass deflagration model for Type la supernovae (SNe la), a persisting free 
parameter is the initial morphology of the flame front, which is linked to the ignition process in the progenitor white dwarf. 
Previous analytical models indicate that the thermal runaway is driven by temperature perturbations (“bubbles”) that develop 
in the white dwarf’s convective core. In order to probe the conditions at ignition (diameters, temperatures and evolutionary 
timescales), we have performed hydrodynamical 2D simulations of buoyant bubbles in white dwarf interiors. Our results show 
that fragmentation occurring during the bubble rise affects the outcome of the bubble evolution. Possible implications for the 
ignition process of SNe la are discussed. 
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1. Introduction 


Advances in computer simulations of type la supernovae 
(SN la) have led to a rough consensus over the nature of 
these explosions, even if several issues remain controver¬ 
sial. SNe la are believed to be the outcome of thermonu¬ 
clear explosions of carbon-oxygen white dwarfs (CO WDs) 
in binary systems. Among po ssible explosion scenarios (see 
iHillebrandt & Niemevetll20()Cll for a review of them) the sin¬ 
gle degenerate scenario fares best in exp laining the observed 
homogeneity of SNe la (lLeibundgull2000h . In this model a CO 
WD approaches the Chandrasekhar mass by mass transfer from 
a non-degenerate companion and is disrupted by a thermonu¬ 
clear explosion. 

In this paper we wish to study the processes that lead 
to the ignition of SNe la. The mechanism that ignites the 
thermonuclear flame in the WD core is the link between the 
late stages of the progenitor evolution and the development 
of the initial flame morphology. From a practical point of 
view, it provides the initial conditions for the explosion sim¬ 
ulations. A large set of different initial conditions has been 
tested in the literature of SN la simulations, with results ranging 
from mildly energetic (|Re.in ecke_.et af 2002^ Gamezo et alJ 
l2003t iTravaglio et alJ 2004t iRo nke & Hillebrandl I200-5I) to 


failed explosions ilCalder et al Jl2fl04li in pure deflagration ex¬ 
plosion models. The initial flame setup in multi-dimensional 
simulations of SN la explosions is considered basically 
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a free parameter, though conshained by some analytical 
works lln arcfa-Senz & W ooslevI 1199.4 IWooslev et alJ 120041: 
IWiinsch &^4rosie^1l20oT Sect. l2.2> . Consequently most sim¬ 
ulation schemes are tested against different sets of initial 
flame shap es, allowing a comparison of the results for yari- 


ous setuDS jReinecke et alJ 

20024 Ronke & Hillebrandtll200.5l; 

iGarcla-Senz & Brayol2005 

iLiyne et al.|200.5j). 


Among the possible initial conditions the most intuitiye 
one, also for it s use in ID simulations, is the centrally ig¬ 
nited flame (e .g.lNiemeyer & Hille brandfll 99.51 iReinecke et alJ 
ll999il2002aHC.ame7o et alJl2003li . A different ignition condi¬ 
tion is proyided by a number of spherical flame kernels placed 
in the innermost p art of the WD and detached from the center . 
Early works in 2D jNiemeyer et all 1996t IReinecke et alJl999l) 
impleme nted only a few blobs per quadrant, while new er sim¬ 
ulations jReinecke et alJl2002bt iTrayaglio et alJE004l) haye a 
finer spatial resolution and allow the implementation of more 
seeds. Neyertheless, the maximum number of blobs that con¬ 
stitute the initial flame is still constrained by resolution rather 
than physics. This situation ma y be cured by the use of co¬ 
expanding computational grids dRopkd 1200-511 . A hybrid ap¬ 
proach between centrally ignited flame and m ulti-spot scenario 
is presented by iRopke & Hillebrandl il2004l in a full-star 3D 
simulation. 

The influence of different initial conditions on the outcome 
of the explosion has been discussed in seyeral papers. There 
is general agreement that an initially larger number of blobs 
produces more yigorous explosions because the initial flame 
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surface is relatively larger and can provide more flame accel- 
eration. Also the compariso n of different initial conditions in 
Ropke & Hillebrandlj feoOSi) shows that the hybrid “foamy” 
initialization provides more seeds for the development of in¬ 
stabilities and gives a larger total energy and production of 
^®Ni with respect to the model ignited at the center. In the 
same work, the com parison with the failed explosion found by 
ICalder et al.l ll2004 is addressed. Their simulation starts from 
a single spherical flame seed, displaced slightly off-center. 

Since the features of SN la explosions depend on the num¬ 
ber of igniting points in the WD core, an interesting question 
pertains to an estimate of this number. In the papers cited above 
the study of this problem was conducted by means of an a pos¬ 
teriori evaluation of the explosion features as a function of the 
initial conditions, without exploring t he physics of the igni tion 
process itself (with the exception of IWooslev et '^ 12004 . A 
successful ignition theory has to provide a physical basis for the 
choice of favored initial conditions in SNe la explosion models. 

In this paper, a new approach to the ignition problem is 
based on 2D simu lations performed with the FLASH code 
jFrvxell et alJ200ol) . We will use an indirect approach; the fea¬ 
tures of the ignition process will be explored by studying tem¬ 
perature perturbations in the WD core (“bubbles”) that are sup¬ 
posed to trigger the ignition (cf. Sect. 12.2> . Using the results of 
the bubble simulations, we are able to obtain important clues 
about the ignition physics. 

The paper is structured as follows; in Sect. |2lthe basic ideas 
of the physics of ignition are introduced. Section 0 describes 
the setup of the simulations. Sectionl^focuses on the bubbles, 
presenting the relevant physics and the issues related to the 2D 
simulations. Some important parameters for the bubble evolu¬ 
tion are highlighted and their role is explored in Sect. |3 In 
Sect. Ihjwe use the results of the simulations in order to study 
possible implications for the ignition of SNe la. Our conclu¬ 
sions are summarized in Sect.0 

2. The ignition process of SNe ia 

2.1. Final stage of the progenitor evolution 

According to stellar evolution models, CO WDs are the end¬ 
points of the e volution of main seq uence stars in the mass range 
M » 3-9Mq jUmeda et al.ll999l) . The energetics in the WD’s 
interior during the accretion phase is determined by the follow¬ 
ing processes; 


mechanism describe d for example bv lClavtonllll984 and 
IWinget et alJ ( |2004 . As a result of accretion, the density 
increases above 2 x 10^ g cm“^, the thermal photons do not 
have enough energy for plasmon excitation and the neutrino 
production in the WD core is strongly suppressed. 

In the e volutionary plane of cen tral density and temperature 
(pc ; Tc) (cf. lYoon & Langeil2003h . the locus where the energy 
release due to carbon burning is equal to the energy loss by 
neutrino emission is traditionally referred to as the “ignition 
line”. Actually it marks the start of the last phase of the pro¬ 
genitor evolution. In order to get rid of the energy output of 
the carbon burning, the WD develops a convective core whose 
evolution is crucial for the knowledge of the ignition proper¬ 
ties. T he duration of the convective phase is of the order of 1 0^ 
years (iHillebrandt & NiemeveTl200r)HYoon & LangeTl2004 . 

When the core temperature reaches 6 - 7 x 10^ K, the con¬ 
vective zone encompasse s most of the WD ma ss. Convection 
enters a reactive regime jWooslev et al.lEoO^) . in which the 
convective turnover timescale and the nuclear timescale Tn, de¬ 
fined as the time required to significantly reduce the carbon 
abundance in an isolated region of the fl uid, are comparable. 
At T ^ 10® K the critical temperature (iTimmes & Woosle^ 
Il992h for carbon burning is reached, i.e. the temperature where 
the energy generation rate is equal to the heat conduction rate. 
This point defines the start of the SN la explosion, the ignition 
of the thermonuclear runaway. 

Beyond the description of this sequence of events, a con¬ 
sistent ignition model of SNe la is required to answer more de¬ 
tailed questions. In particular, the initial location and morphol¬ 
ogy of the flame in the WD, the time evolution of the ignition 
process and the thermodynamic state of the WD core prior to 
runaway need to be explored. One possible approach to study 
the multi-dimensional features of reactive convection is a di¬ 
rect numerical simulation of this phase. A 2D simulation per¬ 
formed with an implicit code and limi ted to the last few hour s 
before the runaway was performed by Hoflich & Steinlll2n02li . 
According to the results of this work, the ignition occurs at the 
center of the WD and is triggered by compressional heating 
in the convective flow . These conclusions were challenged by 
IWooslev et al.l ll2flfl4 who claim that the central ignition was 
either an artifact caused by the forced symmetry of the 2D sim¬ 
ulation, or caused by insufficient resolution. 

2.2. Temperature fluctuations in the WD core 


- Compressio nal heating c aused by accretion, described in 
detail e.g. bv llbenllll982h . 

- Nuclear burning; Compression and heating trigger hy- 
drostatic carbon burning in the WD core. According to 
iNomotol lll982l) the burning starts in the pycnonuclear 
regime and turns to the thermonuclear regime for T > 
5 X 10^ K. 

- Neutrino emission; The energy loss by neutrino emi ssion is 
important at densities below about 2 x 10® g cm~^ jArnetJ 
Il97ll; iNomoto et alJIl9^ IWooslev & Weav^ll986l) . In 
these conditions, the dominant contribution to neutrino 
losses comes from the production of plasmon neutrinos, a 


The idea that the SN la ex plosion is ignited by floating bu bbles 
has been proposed first bv lGarcfa-Senz & WooslevU 19951) . The 
formation of these temperature fluctuations in the WD’s con¬ 
vective core is connected to the turbulent behavior of the WD 
fluid and can be explained in terms of velocity fluctuations. 
IWunsch & Woosl^ ll2004 show that, when a fluid element, 
during the convective motion, approaches the energy generat¬ 
ing core of the WD' with a velocity smaller than the average 
convective velocity, it stays there for a slightly longer time than 

' Considering the h igh sensitivity of the C -buming rate on tempera¬ 
ture, it may be proved Iwooslev et alJ2004 that one half of the energy 
is produced in a small (estimated mass of 0.01 Mq, corresponding to 
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average and becomes hotter (and less dense) than the surround¬ 
ing material. Therefore, the bubble is subject to buoyancy and 
rises from the center of the WD. 

The rising bubbles play a key role in the ignition because 
they are regarded as the “seeds” for the subsequent flame prop¬ 
agation. The initial flame in the progenitor is determined by 
motion of the bubbles inside the convective core. The study of 
the bubble features is an important part of investigating the ig¬ 
nition process. 

iGarcfa-Senz & WooslevI ll 19951) have developed an analytic 
model for the buoyant evolution of burning bubbles. They have 
performed a parameter study in which they investigate the bub¬ 
ble motion as a function of the initial bubble position, the di¬ 
ameter, and the initial temperature excess. The conclusion is 
that the bubbles reach the runaway temperature when they are 
at a central distance in the range 100 - 250 km, moving with a 
speed of about 100 km s'. 

The topic has been revisited by IWooslev et al 1 il2004 . On 
the topic of ignition they conclude that it is initiated by bubbles 
at a central distance of up to about 150 - 200 km, when the 
central WD temperature is about 7.7 - 7.9 x 10® K. Moreover, 
they claim that a multi-point ignition is possible. This results 
from the fact that the e-folding timescale for the number of 
bubbles is comparable with the time for the expansion in the 
SN la explosion to shut off the ignition (~ 0.1 s). 

Finally, the ignition has been addresse d in a more general, 
heuris tic model of turbulent convection bv IWunsch & WooslevI 
Two convective flow patterns are considered, an 
isotropic turbulent flow and a dipolar jet flow. As far as the 
central WD temperature is concerned, the conclusions for the 
SN ignition are similar to previous work. On the location of 
the hottest points in the WD the authors consider the compe¬ 
tition between nuclear heating and adiabatic cooling and esti¬ 
mate that the ignition points should be concentrated at about 
100 km from the WD center. For an isotropic flow it means 
that the ignition points should occur in a shell of 100 km of di¬ 
ameter while for a dipole flow these points would be clustered 
rather along the flow axes. 


3. Numerical tools 


The 2D simulations which we will present here have been pro¬ 
duced using FLASH (version 2 .3), a parallel adaptive-mesh 
hydro code jFrvxell et alJl200oil . The FLASH code has been 
used, among other topics, for the stu dy of SN la explosions 
llCalder et alJl2004l IPlewa et '^120041) and, in a different con- 
text, also in the study of the morphology o f rising bubbles 
jRobinson et alJ2004t Brflggen et al.l2005alhl) in cooling flows 
in galaxy clusters. 

Figure^shows the computational domain at the beginning 
of a 2D bubble simulation. As explained below, the computa¬ 
tional domain encloses only a small part of the WD. In this 
domain the bubble is set at rest. Since most of the calculations 
are performed with bubbles of 1 km in diameter, the following 


about 130 km for the WD model examined by those authors) central 
part of the convective core. 



0 1x10= 2x105 3x105 4x105 5x10= 

X (cm) 


Fig. 1. Temperature plot, that shows a part of the initial setup 
of a bubble simulation. The temperature is coded in gray-scale, 
in Kelvin. The grid is superimposed on the plot; each square is 
defined as a “block” of the AMR structure and contains 8x8 
computational cells. The plot refers to a calculation with an 
initial bubble temperature of 7.7 x 10® K, bubble diameter 1 km 
and initial distance from the center 100 km. 


description will refer to this case (for different bubble diame¬ 
ters one just has to scale the lengths accordingly). 

The size of the computational domain is 5 x 20 km. The 
bubble is initialized as a temperature perturbation in pressure 
equilibrium with the surrounding matter. The values of the ther¬ 
modynamical variables in the domain are obtained from a one¬ 
dimensional WD model provided by S. Woosley, with mass 
Mwd ~ 1.38Mq, radius Rwd ~ 1600 km, Tc - 7 x 10® K, 
Pc = 2.55 X 10® g cm The extent of the convective zone in 
the model is about 1000 km. 

Since the computational domain encloses a relatively small 
part of the WD, a reasonable approach is to neglect the effect of 
curvature and to map the data from the WD model in a plane- 
parallel approximation. The values of physical quantities on the 
computational grid at coordinate y are taken from the quantities 
in the WD model at the radius R-y^+y, where R is the param¬ 
eter which expresses the initial distance of the bubble from the 
WD center and yt = 2.5 km is the y coordinate of the bubble 
center in the computational domain. Before initialization, the 
pressure and density data are slightly modified to ensure the hy¬ 
drostatic equilibrium, as described bv IZingale et alJ ll2002i) . A 
spatially constant gravitational force, computed from the WD 
model, is applied to the computational domain. This assump¬ 
tion is not strictly valid because the variation of the gravita¬ 
tional acceleration along the whole computational domain (20 
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km) is relevant. However, we do not expect this to affect our 
results since the extent of the bubble motion is always smaller 
than the computational domain, namely of the order of a few 
kilometers (Fig. |5}, comparable with the bubble diameter. The 
boundary conditions are set to be reflecting in the upper and 
lower boundaries, and periodic in the x-boundaries. 

The simulations are performed in 2D Cartesian geometry. 
In principle, this is not the most natural choice of geometry 
since a bubble is a spherical object and it would be better rep¬ 
resented in 2D in cylindrical coordinates, exploiting the ax¬ 
ial symmetry. Nevertheless, some tests have been performed, 
both, in 2D Cartesian and cylindrical geometry. Simulations in 
cylindrical geometry present special difficulties, with respect to 
the Cartesian geometry, because of larger numerical dispersion 
(cf. Sect. I4.3> and the appearance of “axis jets” which distort 
the bubble morphology. The tests show that the choice of the 
Cartesian geometry is acceptable for following the thermal evo¬ 
lution of bubbles. 

The adaptive mesh refinement implemented in FLASH is 
apparent in Fig. [Q Resolution tests have shown that a good 
compromise between computational cost and adequate resolu¬ 
tion is to use five levels of refinement yielding an effective grid 
size of [256 x 1024] zones (cf. Sect. 14.3> . This corresponds to 
a spatial resolution of 2 x 10^ cm at the most resolved level. 

Among the different equations of sta te implemented in 
FLAS H, the Helmholtz EOS described by iTimmes & Swest^ 
1I2OOOI) was chosen for our setup. In order to follow the hy¬ 
drostatic carbon burning the small reaction network iso7 was 
used. This a-network is adequate for the physical prob lem un¬ 
der examination, as confirmed bv lTimm^ et al.l( l2000l) . 

To date, in the latest publicly available version of FLASH 
(v. 2.5) there is no treatment of flame propagation. Therefore, 
the simulations are followed until T 10^ K, which is when 
the flame is going to start, and are then stopped. Since this work 
is devoted to the study of the progenitor’s evolution before the 
runaway starts, this limitation does not concern us much. 

4. Two-dimensional simulations of rising bubbles 

4.1. List of the performed simulations 

In order to cover the broad range of bubble parameters that is 
likely to exist inside the WD, we present a study that explores 
the dependence of the bubble evolution on the three main pa¬ 
rameters; 

- The initial bubble diameter, D. 

- The initial bubble temperature, T. 

- The initial distance of the bubble from the WD center, R. 

The aim of the study is to understand under which (phys¬ 
ically meaningful) conditions thermonuclear runaway ensues 
and what the related timescales are. 

The thirteen calculations performed for our parameter 
study are listed in Table [2 In the first column the central dis¬ 
tances of the bubbles are given, together with the background 
temperature of the WD at that distance from the center (in units 
of 10* K). In the other four columns, the simulations are iden¬ 
tified by the initial temperature, again in units of 10* K. The 


Table 1. Scheme of the simulations which have been performed 
for the parameter study. 


Distance from the 

WD center and 
background temperature 



Parameters 

of the simulations 


R = 50 km 

7.3 

7.5 

7.6 


Background = 6.96 








7.7, D = 5 km 


R = 100 km 

7.3 

7.5 

7.7 

7.9 

Background = 6.83 



7.7, D = 0.2 km 


R = 100 km 

7.3 

7.5 

7.7 

7.9 

Background = 6.63 






bubble diameter D is fixed at 1 km unless explicitly specified. 
Simulations that reach the thermal runaway are in boldface. To 
sum up, the initial parameters have been varied in the following 
way: 

1. D (bubble diameter); for the case with T = 7.7 x 10* K, 
R = 100 km, the three values 0.2, 1, 5 km have been tested. 
The lower limit is constrained by computational feasibil¬ 
ity, decreasing the timestep of the simulations together with 
the spatial resolution because of the CFL condition for the 
equations of hydrodynamics. 

2. T (bubble temperature); 7.3 - 7.9 x 10* K. The upper end 
of this range is set in order to avoid too extreme tempera¬ 
ture contrasts with respect to the background. For bubble 
temperatures smaller than the lower end of the range, the 
nuclear timescale is longer than about 5 s. It seems unlikely 
that a bubble can have such a long evolution, without being 
disrupted as described in Sect. 14.31 

3. R (distance from the WD’s center): the three values 50,100, 
150 km have been explored. Following IWunsch & WooslevI 
ll2004ll . we assume that the bubbles are produced in the en¬ 
ergy generating core of the WD, which extends approxi¬ 
mately over this range of central distances. 

4.2. Diagnostic quantities 

First, we will define some quantities that are useful for the in¬ 
terpretation of our simulations. 

A first difficulty in defining these quantities is the lack of a 
good criterion about the bubble location in the computational 
domain. This is of course a typical problem, when dealing with 
an Eulerian system of coordinates. In principle, hot material 
(the bubble is set to be rather hotter than the background) could 
keep track of the bubble evolution. Unfortunately, in Sect. B31 
we show that numerical diffusion hampers this approach and 
one cannot make safe use of “bubble averaged quantities”. 
Hence, the description of the thermal evolution of the bubble 
will be done by studying the maximum temperature in the bub¬ 
ble. 

A useful quantity for the following analysis of the bubble is 
its area. As written above, it is difficult to provide an unambigu- 
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ous definition of this variable. Test simulations have established 
that a good criterion is given by the following expression: 

'l^bubble “ zone ( 1 ) 

bubble 


where the sum is calculated over the zones where Tzone > 
^threshold and the threshold temperature is defined by 


^threshold 


0.95 • Tbubbledni for Tbubbledni > 7.5 X 10* K) 
0.96 • Tbubbledni for Tbubbledni < 7.5 X 10* K) 


where Tbubbiedni is the initial bubble temperature. The quantity 
■^zone is geometry dependent. In Cartesian geometry it can be 
simply identified with the zone area ^zone - dx ■ dy, where dx 
and dy are the zone sizes in the two directions. 

It is also useful to introduce the effecti v e gravitational 
accele ration dOarcfa-Senz & WooslevI Il995t IWooslev et alJ 
I 2 OO 4 I) . The acceleration geft exerted by buoyancy on a rising 
bubble is proportional to its density contrast: 


Ap 

gi:«ir) = g(r) — = 
P 


GM(r) ^ AT 

,2 T ■ 


(3) 


In the previous equation, the density contrast between the 
bubble and the background Ap/p can be expressed in terms of 
the temperature contrast AT /T times (5p, the logarithmic deriva¬ 
tive of density with respect to temperature at constant pressure. 


4.3. The physics of the bubble 

There exists an extensive literature on the fluid mechanics of 
rising bubbles. Nonetheless, we are not aware of any work that 
addresses the peculiarities of the present setup (degenerate mat¬ 
ter, very small density contrasts and nuclear burning). 

At first, our discussion will focus on the following “refer¬ 
ence choice” of parameters: initial temperature T = 7.7x 10* K, 
initial diameter D - \ km, initial central distance /? = 100 km 
(Fig. m These values have been chosen to lie well inside 
the range of variation of the three parameters, dis cussed in 
Sect. 14. H and also suggested bv lWooslev et al.l l l2004l) . First we 
will describe the results for these reference parameters, while 
in Sect. |3we will explore the results of varying these parame¬ 
ters. 

If we neglect adiabatic cooling (as motivated in Sect. 16. H . 
the evolution of the bubble depends only on the interplay be¬ 
tween nuclear burning and hydrodynamic al instabilities. 

Hydrostatic carbon burning sets the timescale for the bub¬ 
ble to reach the thermonu clear runaway. This nuclear timescale 
Tn dWooslev et al.l l2004l) has been introduced in Sect. IQ 
Using the evolutionary timescales of the simulations that reach 
a runaway (cf. Table 0 and adding data from further tests an 
analytical fit to Tn has been found. The result is 



Fig. 2. Series of temperature plots, showing the evolution of a 
bubble in a simulation with the initial parameters T - 7.7 x 
10* K, D = 1 km, = 100 km. The temperature is coded in 
gray-scale, in Kelvin. 


The bubble motion is a special case of a system subject 
to the Rayleigh-Taylor instability (in the following, RTI). In a 
gravitational field with the vector of gravitational acceleration 
g pointing downwards along the y-axis, the hotter and lighter 
fluid accelerates upwards. The development of the RTI and the 
features of the motion in this situation have been addressed 
in several works dP avies ; & TavloHrioiiiiTavloidll950t IT ayzerl 

ll9.5.5Uniimm & I.illl988ill.ill996h . A crucial parameter in this 

problem is the Atwood number Af, 


Af = 


P2 - Pi 
P2 + Pi 


(5) 


where pi and p 2 are the bubble and background densities, re¬ 
spectively. In our bubble problem, the density contrast is very 
small because of the degeneracy of the WD matter, and Af is 
of the order of 10 Theory predicts that the rise velocity of 
the bubble tends to a value Vb that is determined by the equilib¬ 
rium between buoyancy and drag forces. An analytical model 
developed for nonli near, single-mode, classical RTI at arbitrary 
Atwood number bv iGoncharovI d2002ll provides the estimate 


Tn 


^ 10 




S 


(4) 


Vb 


2Af j_ 
1 4- Af Ck 


(6) 


p/(10® gem-*). This fit of 


where Tg = ^/(lO* K) and pg 
Tn is in rough agreement with the derivation bv lWooslev et al 
d2004l) for the conditions of interest in the WD core. 


where k ~ 27r/(D/2) is the wavenumber of a perturbation 
of size of order of the bubble radius Z)/2, and C is a nu¬ 
merical constant whose value is 3 in 2D and 1 in 3D. From 
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Fig. 3. Evolution of the maximum bubble temperature in simu¬ 
lations performed with the reference choice of parameters and 
increasing resolution, indicated by the different levels of refine¬ 
ment. 

Eq. ® one gets, for the reference choice of bubble parame¬ 
ters, Vb - 1-3 X 10^ cm s '. This is an appropriate order-of- 
magnitude estimate for the bubble velocity in the simulations, 
even if the morphology evolution of the bubble prevents a more 
quantitative analysis. 

The vortical motions produced during the bubble rise cause 
its fragmentation as it is clear from Eig. |2 In the nonlinear 
stage of the RTI a large variety of phenomena c an account for 
the ob served fragmentat ion, as described e.g. bylAbarzhi et al.l 
1 I 2 OO 3 I) and reviewed bv iNakai & Takah^ ill996ll in the frame¬ 
work of inertial confinement fusion. We will not address a de¬ 
tailed study about the nature of these secondary instabilities. 
The features of the fragmentation process and some related nu¬ 
merical issues are described below. 

An approximate estimate for the timescale of the fragmen¬ 
tation process can be given by defining a dispersion timescale 

D 

Tdisp = — (7) 

Vb 

which, with the reference values of the parameters and the use 
of eq. i), gives a timescale of about 0.8 s, a sort of “bubble life¬ 
time”, to be compared with the morphological behavior shown 
inEig.|2l 

During the fragmentation, the typical length scale of the 
bubble fragments decreases in time. Erom a physical point of 
view, one can expect that the dispersion goes on until the typ¬ 
ical length scale of a bubble part is comparable with the mini¬ 
mum size /Imin, defined as the diameter of the bubble in which 
the heat generated by nuclear burning is balanced by heat diffu¬ 
sion. An evalu ation of this quantity i s pe rformed, with different 
deriva tions. bv lWooslev et al.l ( l2004l) and lGarcfa-Senz & Bravol 
ll2005l) . and gives results in the range 10 - 100 cm. 

Eor bubbles whose size is larger than Tmin, the nuclear heat¬ 
ing is larger than the heat loss by conduction. In these bubbles 
the temperature increases (in this phase, in principle, they can 



Fig. 4. Evolution of the normalized bubble area for the refer¬ 
ence choice of parameters, in a comparison among simulations 
with different resolution. 


Fig. 5. Series of temperature plots referring to the same simu¬ 
lation as Fig. |2l showing the evolution of a bubble detail. The 
contours in the plots mark the bubble area, according to the 
definition of Eq.[2 

still reach 10® K and trigger the thermonuclear runaway), un¬ 
til disruption causes their typical size to decrease to Tmin. Then 
they begin to cool down. Thus, the trend of the maximum tem¬ 
perature shown in Fig.|3can be easily interpreted; the bubble is 
heated until either runaway occurs or the bubble is disrupted. 

One should note, however, that the length scale Ttnin for dis¬ 
persion driven by heat conduction is much smaller than the spa¬ 
tial resolution of the simulations (2 x 10^ cm). Tmin is so small 
that a direct simulation to that scale, even with AMR, is not 
feasible. What is actually observed in the bubble simulations is 
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not a physical but a numerical dispersion. Figure|3is meant to 
clarify this important point. The plots are zoomed on a small 
detail, a sort of “bridge” linking the top of the bubble with a 
lower vortical structure. This detail becomes thinner and, as its 
width is approximately 2-3 times the spatial resolution, it cools 
off unphysically. It is thus lost from the bubble area according 
to the definition in Sect. 14.21 

The dependence of the bubble fragmentation process on 
resolution is shown in Fig. |3 As expected, the numerical dis¬ 
persion depends on spatial resolution, and the temperature de¬ 
crease occurs earlier for the simulation with coarser level of 
refinement. The resolution which has been chosen for the pa¬ 
rameter study (five levels of refinement) guarantees however to 
follow reliably the fragmentation process, as shown in Fig@] 
Interestingly, the area evolution in the two most resolved sim¬ 
ulations does not show any significant dependence of the bub¬ 
ble area decrease on the level of refinement of the simulation. 
Probably, this happens because the hydrodynamical instabili¬ 
ties in the more resolved calculation produce structures that are 
globally more numerous, though smaller in size. 

As described above, electron conduction is an impor¬ 
tant ingredient in the bubble physics, but can be neglected 
at the length scales which are explored in this study. It is 
kno wn from the theory of th ermonuclear combustion (see 
e.g. Irimmes & Woosl^l 19921) that heat exchange can inhibit 
the growth of instabilities whic h are smaller than a minim um 
length scale Tcond- According to lTimmes & Wooslevl(ll992l) . 

dcond = 47rv^^g-^j (8) 

where g is the modulus of the gravitational acceleration and 
lS.plp is the density contrast of the bubble with respect to 
the background. The flame velocity v can be meaningfully 
replaced, in the case of volume burning, by dmin/Tn- In the 
case under examination (Tg = 7, pg = 2.5) the typical val¬ 
ues of the quantities in Eq. dH) are (cfr. Sect. E2} Tn = 10 S, 
dmin = 100 cm, Ap/p ^ 10“^ and g ^ 10'° cm s“^. This esti¬ 
mate shows that dcond is of the order of 10“^ cm. It is smaller 
than dmin, therefore the described bubble physics is not affected 
by the stabilization. 

4.4. Outcome of the bubble evolution 

The following overview sums up the main features of the bub¬ 
ble evolution: 

- The initial bubble temperature determines the nuclear 
timescale (Eq. |^. This is the upper limit for the dura¬ 
tion of the bubble evolution. 

- The bubble is fragmented during its motion on a timescale 
approximately given by r^isp (Eq.0l. Physically, the disper¬ 
sion proceeds down to length scales where heat conduction 
is effective in dissipating the energy generated by nuclear 
burning. 

- The previous argument is made more complicated by nu¬ 
merical diffusion. It is not possible to quantify exactly the 
whole duration of the physical dispersion phase described 


previously because this process can be followed in the sim¬ 
ulations only until the typical length scale of the bubble 
substructures are comparable to 2-3 times the spatial reso¬ 
lution. 

The competition between nuclear heating and dispersion is 
the key to understand the outcome of the bubble evolution and, 
consequently, also to the ignition process of SNe la. Roughly 
speaking, if Tn is smaller than Tdisp the burning prevails and the 
bubble reaches the thermonuclear runaway. Conversely, if Tn > 
Tdisp the bubble is disrupted and cools down. Analytical studies 
on ignition that is driv en by floating bubbles have not explored 
the bubble dispers ion dOarcfa-Senz & Wooslevll995[ [Wooslevl 

l2nniHWooslev et al 12004 IWiinsch & Wooslevi2nn4l) . 

4.5. Background turbulence 

In all our simulations the background state of the WD was as¬ 
sumed to be “quiet” and in hydrostatic equilibrium. However, 
the convective flow in the WD prior to runaw ay is turbulent, 
with a n estimated Reynolds number about lO''* dWooslev et al.l 
I 2 OO 4 I) . The integral length scale L, at which the turbulent ki¬ 
netic energy is injected, is given by the press ure scale heigh t, 
about 450 km for the WD model used here dWooslevIbOOll) . 
Assuming Kolmogorov scaling for the typical velocity of the 
convective eddie of length scale I, the characteristic turbu¬ 
lent velocity U at this length scale is expressed by t/(/) » 
U{L){llLyi\ 

In our case, / = 2 x 10^ cm (spatial resolution of the sim¬ 
ulation at the most resolved level) and UiL) is in the range 
50 - 100 km s ', from estimates based on mixing length the¬ 
ory. Therefore, the typical turbulent velocity U{1) is in the range 
1.7 - 3.5 km s ', which is of the same order of magnitude as 
the typical rise velocity of the bubble, according to Eq. (|^. A 
test simulation, performed by imposing a random velocity field 
with an amplitude of the order of t/(/), shows that this stir¬ 
ring affects the bubble evolution only marginally, contributing 
to a moderately enhanced bubble disruption. We conclude that 
this turbulence-induced bubble disruption can be considered a 
second-order effect in the description of the relevant bubble 
physics discussed in Sect. 10 

Our simulations are affected by numerical dispersion and 
hence overestimate the bubble disruption and the subsequent 
bubble cooling. Erom a physical point of view one expects that 
the bubble temperature would increase for a longer time than 
is shown, for example, in Eig. |3 However, the numerical dis¬ 
persion may not be too unphysical because it may mimic the 
effect of the turbulent-induced bubble disruption. A qualitative 
comparison of these two effects cannot be performed with our 
computational tools. 

5. The parameter study 

5.1. The effect of the bubble diameter D 

As outlined in Table [I] the role of different bubble diameters 
was studied on the basis of three calculations, with initial pa¬ 
rameters T - 7.7 X 10^ K and R — 100 km. The initial diameter 
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time [s] 


Fig. 6. Comparison of the evolution of maximum temperature 
in simulations with R - 100 km, T - 1.1 x 10® K and D 
indicated in the legends. 


Fig. 8. Comparison of temperature evolution in simulations 
with = 100 km, D - \ km for four different initial tem¬ 
peratures. 



Fig. 7. Same as Fig.|6l but the evolution of the normalized bub¬ 
ble area is shown. 

D was varied by a factor of five above and below the reference 
value, i.e. D - 0.2,1 and 5 km. The extent of the computational 
domain and the spatial resolution were scaled accordingly in 
order to resolve the initial bubble by an identical number of 
zones. Figures 13 and 0 present a comparison of the maximum 
temperature and the normalized bubble area, respectively. 

As expected, the bubble diameter affects the dispersion 
timescales. Indeed, the dispersion timescale (Eq. 0 scales as 
D F2 since, from Eq. •EJi, Vb oc D F2 as well. This scaling is 
nicely confirmed by the comparison of the area evolution in 
Fig.El While in the simulation with D - \ km the area starts 
to decrease at approximately t ^ 0.8 s, one can see that for 
the bubble with D - 0.2 km the analogous decrease starts at 
f 0.8-(0.2/1)F2:^ 0.35 s. 

Eor the same reason, the bubble with D - 5 km has a larger 
Tdisp. Since the nuclear timescale does not depend on the bub¬ 
ble diameter, this leads to Tdisp > Tn, and the bubbles goes to 
runaway. 



Fig. 9. Comparison of the normalized bubble area in the param¬ 
eter study on temperature. The initial bubble temperatures are 
indicated in the legend, in units of 10® K. 


From this study, one can conclude that t he larger bub¬ 
bles ar e more likely to go to runaway. However. lWooslev et alJ 
il2004 find that turbulent dispersion limits the bubble size in 
the WD to Z) ~ 1 km. 

5.2. The effect of the bubble temperature T 

The initial temperature T is the most interesting parameter be¬ 
cause its role in the bubble physics is manifold. As shown in 
table [3 several values have been tested. For sake of simplic¬ 
ity, a discussion will only be presented for the simulations with 
R - 100 km, D - I km and T ranging from 7.3 to 7.9 x 10® K, 
keeping in mind that the inferi'ed trends are valid also for cal¬ 
culations with other central distances. Figures|3and|3present a 
comparative analysis of the evolution of maximum temperature 
and area. 
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Fig. 10. Comparison of the evolution of maximum temperature 
in simulations with T = 7.5 x 10^ K, D = 1 km and R indicated 
in the legends. 



Fig. 11. Same as Fig. cni but the evolution of the normalized 
bubble area is shown. 


The nuclear timescale depends on the bubble temperature 
in a rather steep way (Eq. 1^. Also the dispersion timescale 
Tdisp depends implicitly on T, approximately via the square 
root of the Atwood number At (cf. Eq. |^. At varies in the ex¬ 
plored temperature range from 2.2 x 10“"^ (T - 7.3 x 10* K) to 
5.3 X 10“"^ (T - 7.9 X 10* K). In simpler terms, the temperature 
contrast between the bubble and the surrounding material is 
linked to the density contrast and hence to the effective gravita¬ 
tional acceleration (Eq.|3- With equal background temperature 
(i.e. equal central distance R) hotter bubbles experience larger 
accelerations and thus faster dispersion. On the other hand, the 
nuclear timescale decreases even faster with temperature. This 
indicates the existence of a threshold temperature above which 
a bubble goes into runaway. At = 100 km the threshold is 
found to be approximately T = 7.9 x 10* K. 


5.3. The effect of the central distance R 

While the parameter T was shown to be the most interesting 
one for the bubble physics, the central distance R is probably 
the most relevant one for the ignition scenario of SNe la. The 
analysis was performed by comparing three simulations with 
the initial parameters T = 7.5 x 10* K, D = 1 km and R equal 
to 50, 100 and 150 km. 

The increase of R in the calculations affects the effective 
gravitational acceleration for two reasons. Errst, the modulus of 
the acceleration g{R) increases with R in the considered range 
of distances. Second, the temperature prohle of the WD is de¬ 
creasing (the background temperatures at various R are given 
in Table 0. So at equal bubble temperature, the temperature 
(and density) contrast is increasing with R. This implies that 
the effective gravitational acceleration increases outwards, too. 

This effect is evident in the comparison between the 
timescales of the two extreme cases, R = 50 and 150 km, from 
the analysis of Eigs.lTfilandll II The simulation with R = 50 km 
goes to thermonuclear runaway. Its area evolution shows an ex¬ 
pansion (partly due to a problem of the definition of “bubble 
area” when the temperature contrast between the bubble and 
the background is relatively small). On the other hand, as pre¬ 
dicted the dispersion timescales of the other simulations are 
noticeably shorter as R increases. 

The role of an increasing geft is also relevant for the value 
of the threshold temperature at which the bubbles at different 
initial distances from the center go into runaway. As shown in 
tabled the bubble at T = 7.5 x 10* K reaches the ignition 
temperature during the simulation with R - 50 km, but with 
= 150 km even the bubble at T = 7.9 x 10* K is dispersed 
before reaching the runaway. 

6. Discussion 

6.1. Further remarks on the physics of rising burning 
bubbles 

One of the most striking differences between the simulations 
presented here and the hrst a nalyti c study on bubble evolution 
by lOarcfa-Senz & WooslevI d 199.51) is that the latter does not 
take into account the dispersion of the bubbles. Consequently, 
in their work the bubble velocity and the path covered during 
the rise are much larger than in the present study. Moreover, 
without dispersion the rise times are rather long, up to about 25 
seconds. 

The rise velocity of the bubbles in the simulations is of the 
order of 1 km s * (cf. Eq. |^. Because the timescale until the 
bubbles either go into runaway or are disrupted (1-2 s) is so 
short, the bubble evolution occurs mostly near their point of 
origin (Eig. |2j- However, the previous argument neglects the 
convective motions in the WD core. 

Placing a bubble at some distance R from the WD’s cen¬ 
ter and letting it rise is a gross simplification because the 
convective velocity is always much larger (in the range 50 - 
100 km s ', from estimates based on the mixing length the¬ 
ory) than the simulated buoyant velocity. So it is the convec¬ 
tion which determines the motion of the bubbles, rather than 
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the rise velocity. Clearly, the convective flow pattern crucially 
affects the spatial distribution of bubbles at runaway. 

As the bubble is advected by the convection, it may 
travel over a distance that is long compared to the esti- 
mated central distance at ignition, undergoing some cooling 
dWunsch & Wooslevi2004l) . The bubble simulations cannot di¬ 
rectly address this physical process. In our setup, the bubble 
is initially at rest and no convective velocities are imposed on 
the background. As a result the bubble travels only a distance 
that is small compared to its diameter. One can prove that adi¬ 
abatic cooling is negligible in this case, as confirmed by test 
simulations in which the nuclear burning was switched off. 


6.2. Implications for the theory of ignition of SNe la 


In this section, our results are applied to the ignition process in 
SNe la. 

First, a rough estimate of the extent of the size of the ig¬ 
nition zone can be made on the basis of bubble lifetimes and 
convective velocities. Taking the values of 1 s and 70 km s ', 
respectively, one finds that the extent of the motion of the 
bubbles advected in the convective flow is about 70 km. 
Considering that the size of the energy-producing core is of 
the order of 100 km, an estimate for the distance covered 
by a rising bubble before it triggers the ignition is of the 
order of 150 km. The uncertainties in the convective veloc¬ 
ity (in the range 50 - 100 km s“') and in the initial bub¬ 
ble location inside the WD core lead to some dispersion 
around this value, which is of the order of 100 km (a frac¬ 
tion of the pres sure scale height). We note that this estimate is 
consis tent withlGarcfa-Senz & WooslevUlOO'^ . IWooslev et al.l 
I and lWiinsch & Wooslev ( 20041) . 


It is clear that our approach for the study of ignition con¬ 
ditions does not allow any conclusions about the departures 
from central symmetry. This issue has to be studied with other 
numerical tools. A promising way seems to be t he use of 3D 
anelastic simulations, adopted by Kiihlen eT all (l2003l) in the 
study of convective flows in massive stars. Interesting results 
of this technique , appli ed to SN la progenitors, are presented 
by iKuhlen et alJ (l2005h . In any case, the present work high¬ 
lights the importance of the progenitor evolution for the initial 
conditions of the SN la explosion. The convective pattern plays 
a crucial role in the ignition process. The initial flame location 
depends mostly on the position of the bubbles, which, in turn, 
are advected by the convective motions inside the WD. 

Another interesting issue is the initial number of igniting 
points. Though our approach for the study of the ignition pro¬ 
cess is indirect, some useful hints can be obtained. 

In the present parameter study, several values of the ini¬ 
tial temperature have been tested without making assumptions 
concerning the probabili ty distribution functio n (PDF) of the 
temperature fluctuations. IWooslev et al.l(l2004l) investigate this 
problem and indicate two possible PDFs, depending on de¬ 
tails on the convective mixing in the WD core. In their work 
the PDF is either exponential or Gaussian. Without making a 
choice among these two models, from the generation process 
of the temperature perturbations it is intuitive that bubbles with 


a relatively large temperature contrast with respect to the back¬ 
ground are less likely to be generated. 

Here a thought experiment of two different ignition scenar¬ 
ios is useful. In the first case, let a hot bubble (T > 7.5 x 10^ K) 
be in the core of a WD with central temperature - 7.0 x 
10^ K. The estimated nuclear timescale is about 2 s. If this 
bubble is on the “hot tail” of the temperature PDF, one can as¬ 
sume that the probability of generating other bubbles with this 
T (or hotter) within Tn is not large, while colder bubbles would 
not have time to reach the runaway temperature. In this sce¬ 
nario, the ignition would occur in one (or a few) igniting points. 
However, considering that the dispersion in more effective for 
bubbles of a large temperature contrast with the background, 
such kind of ignition scenario is likely to fail. If the explosion 
is not initiated, the WD goes on increasing its central tempera¬ 
ture. 

The other scenario is directly connected with the previous 
idea. In this second idealized setup, let the WD have a larger 
central temperature, > 7.5 x 10* K, and the temperature 
fluctuations be relatively mild with respect to the background. 
Under these assumptions, the estimated Tn of the bubbles is 
about 1 s, a value compatible with multi-spot ignition, since 
the e-folding time of the number of bubbles in the WD core 
is about 0.1s jWooslev et al .12004 . Moreover, these perturba¬ 
tions are not supposed to be very hot compared to the back¬ 
ground. The probability for such bubbles to be generated, ac¬ 
cording to whatever PDF, should be large, allowing the pres¬ 
ence of very many of them in the WD core. 

From these arguments the multi-point scenario emerges as 
a good candidate for an ignition model of SNe la. 

7. Conclusions 

We have studied the ignition process in SNe la by simulating 
the evolution of buoyant bubbles in the WD core. After dis¬ 
cussing the bubble physics, the role of the relevant parameters 
was studied by means of a grid of thirteen simulations. This 
parameter study shows that floating bubbles may either ignite 
almost in place or be dispersed by fragmentation. The previ¬ 
ous results found an application in the study of ignition process 
in SNe la. The discussion in Sect. l6.2l r)rovides some clue that 
the multi-spot ignition model is the favored one. This bubble 
distribution is consistent with the initial conditions assumed in 
many recent 3D simulation s of SNe la explosions. 

As pointed out e.g. by iHillebrandt & Niemeveil ll2000l) . a 
successful model for SNe la should explain, both, the homo¬ 
geneity of the class and, hopefully as a function of one or few 
parameters, the diversity of observational features. The present 
study indicates that the ignition process contributes to the ho¬ 
mogeneity of SNe la. The explanation comes from examining 
the second scenario described in Sect. 16.21 According to the 
PDFs (no matter which one), there is some probability for the 
occurrence of one (or a few) relatively hot bubble among the 
the mild temperature fluctuations. In principle, a very hot bub¬ 
ble could lead to a runaway before the colder ones, resulting 
in premature ignition in one or very few points. The physics of 
the bubbles leads to a sort of self-regulation because the effi¬ 
ciency of the disruption increases with the bubble temperature 
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(Sect. 15. and the hottest bubbles are disrupted more effec¬ 
tively. Therefore, one can conclude that a large range of pos¬ 
sible ignition condi tions is unlik e ly. Th e singly-ignited initial 
m odel proposed by Calder et al.l 1120041) and further analyzed 


bv IPlewa et alJ ll2004 ) could be interpreted as an ignition model 
coming from a single-bubble runaway. This kind of ignition is 
not favored by the previous probability arguments. 

Of course, this homogeneity will only apply if the under¬ 
lying convective patterns in the WDs are homogeneous. The 
diversities in the convective flow dWunsch & Woosl^l2004l) 
are potentially able to affect noticeably the ignition conditions. 
Their role in producing the observed range of diversity in SNe 
la has not yet been explored. 

As far as WD central temperature, bubble temperature 
and radius of the ignition zone are concerned, the results 
for the ignition are essentia lly i n agreement with the flnd - 
ings of lWooslev et al.l (l2004l) andlWiinsph Woosleyl j| 2004l) . 
While our study an d the w ork by Wooslev et al. i2004 ) and 
IWunsch & Woosl^ ) l2004l) have the same theoretical back¬ 
ground, our approach and methods are widely different. 

We discussed the simplified assumptions that this work is 
based on, and their role in the analysis of the results. In partic¬ 
ular, possible limitations are introduced by the uniform back¬ 
ground hypothesis (Sect. 14.51 and by neglecting adiabatic cool¬ 
ing (Sect. 16.l1 . 

In SNe la simulations that implement multi-point igni¬ 
tion, the diameter of the flame seeds (and consequently their 
number) is set by the spati al resolution of the simulation. In 


Ropke & Hillebrandt] 120051) this diameter is 7.0 km, and future 


simulations, performed with better computational resources, 
will be able to resolve progressively smaller scales and allocate 
more bubbles. The dependence of the explosion features on 
the number of bubbles, when this number is larger than about 
100, opens new ins ight on the explosion mechanism of SNe la 
jRopke et al .1120051) . At this stage the increased spatial resolu¬ 
tion will not help to improve the results unless there is a better 
understanding of the link between the progenitor evolution and 
the early explosion phase. Future SN simulations should also 
take into account the short (~ 0.1 s) temporal evolution of the 
ignition process, and its interplay with the ongoing explosion. 
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